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1 Introduction 

Most models of low-scale gravity allow the production of non-perturbative gravitational 
states, such as micro black holes and string balls (highly excited string states) at Large 
Hadron Collider (LHC) collision energies [1~4]. This is due to the fundamental gravita¬ 
tional scale being comparable to the electroweak scale (Mew) in these gravity models. If 
black holes or string balls are produced at the LHC with masses much higher than this 
fundamental gravitational scale, they behave as classical thermal states and decay to a 
relatively large number of high-transverse-momentum (high-px) particles. One of the pre¬ 
dictions of these models is the expectation that particles are emitted from black holes 
primarily according to the number of Standard Model (SM) degrees of freedom (number of 
charge, spin, flavour, and colour states). 

To identify high-px) high-multiplicity hnal states resulting from high-mass objects, a 
suitable variable is the scalar sum of the px of the jets in the event, Rx- A low-Rx control 
region is dehned where the background is expected to dominate over any possible new 
physics signal. A ht-based technique is used to extrapolate from the control region to a 
high-Mx signal region to estimate the amount of SM background. 

This paper is organised as follows. The phenomenology of low-scale gravity relevant 
to the search is briefly described in section 2. In section 3, the main components of the 
ATLAS detector are summarised. The Monte Carlo (MC) simulated samples used for 
the analysis are presented in section 4. In section 5, the trigger and event selection are 
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described. The characterisation of the data and the method nsed in the search are given 
in section 6. Section 7 describes the systematic nncertainties, and the resnlting limits are 
given in section 8. Finally, conclnsions are stated in section 9. 

2 Theoretical background and previous results 

Understanding qnantum gravity is one of the main challenges of modern physics. The hier¬ 
archy problem (the relative weakness of gravity compared to the electroweak interaction) 
may be key to that understanding. Two main paradigms for models involving extra dimen¬ 
sions have been formulated: the Arkani-Hamed, Dimopoulos, Dvali (ADD) proposal [1, 2] 
involving large extra dimensions; and a five-dimensional model with a single highly warped 
anti-de Sitter space [3, 4]. These models have our (3+l)-dimensional world residing on 
a brane, which is embedded in a (4+n)-dimensional bulk with n extra dimensions. The 
effective strength of the gravitational interaction inside the brane is weakened by the large 
volume of the extra dimensions or red-shifted by the warp factor along the extra dimension. 
This weakening of the gravitational strength results in a diminished effective Planck scale 
Md \n the (4+n)-dimensional world, relative to the familiar Planck scale M-p\. In the ADD 
model, there are a number n > 1 additional flat extra dimensions, and is determined 
by the volume and shape of the extra dimensions. 

If M^) ~ MeW) several low-scale gravitational signatures may be probed in collider 
physics experiments. Some of the most interesting are the possible existence of non- 
perturbative gravitational states such as black holes [1-4], string balls [5] (in the context of 
weakly coupling string theory), and higher-dimensional branes. 

Within the context of the ADD model, experimental lower limits on the value of M/j [6] 
were obtained from experiments at LEP and the Tevatron [7, 8], as well as at ATLAS [9] and 
CMS [10], by searching for the production of the heavy Kaluza-Klein gravitons associated 
with the extra dimensions. The most stringent limits come from the LHC analyses [9, 10] 
that search for non-interacting gravitons recoiling against a single jet, and range from 
Mo > 3.1 TeV, for n = 6, toME)>5.2 TeV, for n = 2. Several searches for black holes 
and string balls are also performed by ATLAS [11-14] and CMS [15-17]. 

In proton-proton collisions with centre-of-mass energy ^/s, classical black holes form 
when the impact parameter between two colliding partons, with centre-of-mass energy y/§, 
is less than twice the gravitational radius rg of a black hole of mass equal to y/§ [18, 19]. 
Black holes are assumed to be produced over a continuous range of masses above a certain 
threshold Mth > up to ^/s. Semi-classical approximations used in the modelling are 
valid for masses only well above Mo, motivating the use of a minimal threshold Mth to 
remove contributions where the modelling is not reliable. 

Most low-scale gravity models assume classical general relativity to predict the pro¬ 
duction cross section for black holes {a ~ string balls, and use semi-classical 

Hawking evaporation (a completely thermal process due to quantum effects) to describe 
their decay [20]. The decay process is described by black-body radiation at the Hawking 
temperature (Hagedorn temperature for string balls) with the expectation that the radiated 
particle species are produced according to the number of SM degrees of freedom and are not 
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affected by the strengths of the SM forces. The emissions are modihed by spin-dependent 
qnantnm statistics given by the Fermi-Dirac or Bose-Einstein distribntions. In addition, 
the emissions are modihed by gravitational transmission factors [20] (gray-body factors), 
which depend on the spin of the emitted particle, as well as the angnlar momentum of 
the black hole, and can be sizeable for vector particle emission from rotating black holes. 
Once black holes are produced, they evaporate causing their mass to be reduced with each 
emitted particle. In the context of weakly coupled string theory, black holes transition to 
string balls at a minimum black hole mass Mmin ~ where Mg is the string scale and 

Qg is the string coupling constant [5, 21]. When the black hole mass is reduced to approxim¬ 
ately Md (or Mg for string balls), the black hole is said to be in a remnant state, which is 
expected to only be describable by a theory of quantum gravity. This study only considers 
unstable black hole remnants, and black holes and string balls that are short lived. 

The production and decay of black holes and string balls lead to hnal states distin¬ 
guished by a high multiplicity of high-pT particles, consisting mostly of jets arising from 
quark and gluon emission. Since black hole decay is considered to be a stochastic process, a 
different number of particles, and thus jets, can be emitted from black holes with identical 
kinematics. 

3 ATLAS detector 

The ATLAS experiment [22] is a multi-purpose particle physics detector with a forward- 
backward symmetric cylindrical geometry and nearly dvr coverage in solid angle. ^ The layout 
of the detector is dominated by four superconducting magnet systems, which comprise a 
thin solenoid surrounding inner tracking detectors and three large toroids, each consisting 
of eight coils. The inner detector consists of a silicon pixel detector, a silicon microstrip 
detector, and a transition radiation tracker, with a combined coverage up to \r]\ = 2.5. 
In the pseudorapidity region |r/| < 3.2, liquid-argon (LAr) electromagnetic (EM) sampling 
calorimeters are used. An iron/scintillator tile calorimeter provides hadronic coverage over 
\r]\ < 1.7. The end-cap and forward regions, spanning 1.5 < |p| < 4.9, are instrumented with 
LAr calorimetry for EM and hadronic measurements. The muon spectrometer surrounds 
these, and comprises a system of precision tracking and trigger chambers. A three-level 
trigger system is used to select interesting events [23]. The Level-1 trigger is implemented 
in hardware and uses a subset of detector information to reduce the event rate to at most 
75 kHz. This is followed by two software-based trigger levels which together reduce the 
event rate to about 300 Hz. 


^The ATLAS detector uses a right-handed coordinate system with its origin at the nominal interaction 
point in the centre of the detector and the z-axis along the beam direction. The a;-axis points toward 
the centre of the LHC ring, and the y-axis points upward. Cylindrical coordinates (r, (j)) are used in the 
transverse plane, (j) being the azimuthal angle around the beam pipe. The pseudorapidity rj is dehned in 
terms of the polar angle 6 hy g = — ln[tan(0/2)]. 
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4 Monte Carlo simulation 


All background estimates in this analysis are derived from data. However, SM MC simulated 
events are used to estimate the relative background contributions from different processes 
expected in the data sample, and to develop and validate the analysis methods. 

The dominant background in the search region consists of QCD multi-jet events, with 
small contributions from top quark pair production (ft), 7 +jets, VH+jets, and Z+jets. 
Single-top-quark and diboson processes contribute negligibly to the selected samples. The 
baseline samples of inclusive jets are generated using PYTHIA 8.160 [24] implementing LO 
perturbative QCD matrix elements for 2 — 7-2 processes and px-ordered parton showers 
calculated in a leading-logarithmic approximation. The ATLAS AU2 set of MC parameters 
(tune) [25] and the CTIO [26] PDFs are used with these samples. Herwig+-i- 2.6.3 [27] dijet 
samples with the ATLAS EE3 tune and CTEQ6L1 [28] PDFs, and ALPGEN 2.14 [29] multi¬ 
jet samples hadronised with PYTHIA 6.427 with the ATLAS Perugia 2001C tune and the 
CTEQ6L1 PDFs are used for comparisons. The tt, 7 +jets, VF+jets, and Z+jets samples 
are generated using SHERPA 1.4.0 [30] with CTIO PDFs. All MC simulated background 
samples (except ALPGEN) are using the full GEANT4 [31] simulation. 

Signal acceptances are determined using MC simulated events. Signal samples are gen¬ 
erated using the MC event generators CHARYBDIS2 1.0.2 [32] and BlackMax 2.02.0 [33]. Two 
generators are used since they model the remnant decay slightly differently and neither im¬ 
plements all the models considered in this analysis. CHARYBDIS2 is used to produce samples 
for non-rotating, rotating, and low-multiplicity remnant black holes, and for an initial- 
state graviton radiation model. BlackMax is used to produce samples for non-rotating 
and rotating black holes, and for final-state graviton emission and initial-state photon ra¬ 
diation models. The initial-state radiation is modelled to occur after, rather than before, 
black hole formation. In addition, CHARYBDIS2 is used to produce non-rotating and rotating 
string ball samples. Both generators use a leading-order parton distribution function (PDF) 
MSTW2008 [34], the ATLAS AU2 tune, and the PYTHIA 8.165 generator for fragmentation. 
The most important parameters that have significant effects on black hole production are 
Mth, Mu (Ms for string balls), and n. Signal samples are produced for many values of these 
parameters. The MC simulated signal samples are passed through a fast simulation of the 
ATLAS detector [35]. The fast simulation uses a parameterised response of the calorimet¬ 
ers, and GEANT4 for the other parts of the ATLAS detector. The difference in signal yield 
with respect to a full GEANT4 simulation of the ATLAS detector [36] is negligible. 

Additional proton-proton collisions are modelled by overlaying minimum bias events 
on the simulated signal and background events according to the luminosity profile of the 
recorded data. The MC simulated events are reconstructed and analysed with the same 
procedures as used on data. 

5 Trigger and data selection 

The data used in this analysis were recorded in 2012, with the LHC operating at a centre- 
of-mass energy of ^/s = 8 TeV. All detector elements are required to be fully operational. 
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and a total integrated luminosity of 20.3 fb ^ is used in this analysis with a luminosity 
uncertainty of 2.8%. It is derived following the same methodology as that detailed in 
ref. [37]. 

The events used in this search are selected using a high-T^x trigger, which requires at 
least one jet of hadrons with px > 170 GeV and a high scalar sum of transverse momentum 
of all the jets in the event. The trigger is fully efficient if the event has > 1.2 TeV, as 
required in this analysis. 

Events are required to have a primary vertex with at least two associated tracks with 
px above 400 MeV. The primary vertex assigned to the hard scattering collision is the 
one with the highest XltrackT'T’ "w^Gre the scalar sum of track is taken over all tracks 
associated with that vertex. 

Since black holes and string balls are expected to decay predominantly to quarks and 
gluons, the search is simplified by considering only jets. The analysis uses jets of hadrons, as 
well as misidentified jets from photons, electrons, and r leptons. The incorrect calibration 
of photons, electrons, and r leptons using the hadronic energy calibration leads to small 
energy shifts for these particles, but since a particle of this type is expected to occur in less 
than 0.6% (as determined from simulation studies) of the events in the data sample, they 
do not contribute significantly to the resolution of global quantities. 

The anti-fej algorithm [38] is used for jet finding, with a radius parameter R = 0.4. The 
inputs to the jet reconstruction are three-dimensional topo-clusters [39]. This method first 
clusters together topologically connected calorimeter cells and then classifies these clusters 
as either electromagnetic or hadronic. The classification uses a local cluster weighting cal¬ 
ibration scheme based on cell-energy density and longitudinal depth within the calorimeter. 
Based on this classification, energy corrections described in ref. [40] are applied. Further¬ 
more, jets are corrected for pile-up. The jets are required to have px > 50 GeV and \r]\ < 2.8 
in this analysis. 

6 Background estimation method 

Events are selected if they pass the high-i7x trigger and have i7x >1-5 TeV. The discrim¬ 
inating variable chosen for this analysis is 77x- Figure 1 shows the 77x distributions for 
different inclusive jet multiplicities. Data as well as MG simulations of the most significant 
SM contributions to the 77x distributions are shown. The different MG contributions are 
first weighted according to their cross sections, and then the total SM contribution is nor¬ 
malised to the number of data events in the region 1.5 < i7x <2.9 TeV for each inclusive 
jet multiplicity. As can be seen in figure 1, the expected SM background is dominated 
by QGD jet production. The rest of the background processes contribute less than 3% to 
the total background, and therefore the other contributions are neglected in what follows. 
Figure 1 also shows good agreement between SM expectations and data, which is quantified 
in section 8. 

Each event is characterised by the number of jets Ajet and by the value of 77x- The 
(Ajet, 77x)-variable space is divided into two exclusive Hj' regions for the search, which are 
each further divided by inclusive jet multiplicity Ajet- The two exclusive i7x regions are 
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Figure 1: Distributions of the scalar sum of the px of all jets in the event, for 
different inclusive jet multiplicities iVjet for 20.3fb“^ of collision data and MC simulations 
of SM processes. The uncertainties on the data and ratio points are due to the statistical 
uncertainty of the data only. The SM contributions to the background are normalised 
relative to their nominal cross sections and then the total background is normalised to the 
number of data events in the region 1.5 < < 2.9 TeV for each inclusive jet multiplicity. 
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defined as a control region (1.5 < H^i < 2.9 TeV) and as a signal region (//t > 3.0 TeV), for 
each inclusive jet multiplicity. The control region is utilised to ht the distribution, as no 
resonances or threshold enhancements above SM processes were observed in this region. The 
signal region is the kinematic region in which data are compared to the extrapolation from 
the control region to search for enhancements. The search is divided into six overlapping 
regions of inclusive jet multiplicity: iVjet > 3 to A^jet > 8. The region with less than three 
jets is excluded because non-perturbative gravitational states are unlikely to decay to one 
or two jets [19] and the SM background would be larger with respect to the signal than 
in the higher jet-multiplicity regions. Other ATLAS searches [41] have set limits on this 
low-multiplicity region. 

The SM background in each signal region is estimated by htting a function to the data 
in the corresponding control region and then extrapolating the resulting function to the 
signal region. To ht the LTt distribution, a three-parameter Po,Pi,P 2 empirical function 

dN _ po{l-xy^ 

di^T ^ ^ 

where x = Ht'/^/s, is used. 

The selection of the boundaries of the control region is based on 1) stability of the 
extrapolation into the signal region of the function ht to data with respect to small changes 
in the choice of control region, and 2) minimisation of possible black hole (or string ball) 
signal contamination in the control region. 

The effect of possible signal contamination in the control region was studied. A string 
ball sample (n = 6, Mth = 4.5 TeV, Mg = 1.0 TeV, and gg = 0.4) was used in this study 
since it had the largest fraction of events in the control region. The number of string 
ball events was scaled down (by a factor of 27) until its contribution to the signal region 
was at the current level of detectability (three standard deviations above the expected 
background). It was determined that this level of signal would not affect the ht and its 
extrapolation by more than the statistical uncertainty if the upper boundary on the control 
region is below 2.9 TeV. The value of the cross section used in this study has already been 
ruled out at the 95% conhdence level (CL) [11-17]. 

Although the background estimate only relies on data, the validity of the assumption 
that the ht in the control region can be used to estimate the background in the signal 
region was tested using PYTHIA 8, Herwig++, and ALPGEN MC simulated events. Since the 
multi-jet SM simulated events do not contain a signal, the results of the ht to the entire 
H'j' distribution can be compared to the results of the ht to only the control region and 
extrapolating into the signal region. As an example, the results for PYTHIA 8 dijet simulated 
events is shown in hgure 2. These studies show that the ht extrapolation approximates the 
MC simulated events in the signal region to within 20%. This difference is covered by the 
uncertainties, which are described in the next section. 
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Figure 2: The -ffx distributions, showing a comparison of the full range (FR) fit from 
= 1.5 TeV to the last predicted data value with the extrapolation of the fit from the 
control region (CR) 1.5 < H'l < 2.9 TeV into the signal region (i/x > 3.0 TeV) for PYTHIA 8 
simulated events. The uncertainty band includes all the uncertainties described in section 7. 
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7 Systematic uncertainties 


In addition to the statistical uncertainties from the limited number of events in the control 
region, systematic uncertainties arising from the choice of control region and the choice of 
fit function are considered. To estimate the effect of limited number of events in the H^' 
distributions on the fit, pseudo i7x-distributions with the same number of events as data 
are generated using the fit to data in the control region as a probability density function. 
Each pseudo i^x-distribution is fit, and the predicted number of events in each TTt bin is 
calculated and subtracted from the prediction from the fit to data. A distribution of these 
differences between the fit to data and the fit to pseudo distributions is used to derive an 
uncertainty on the fit in each Tfx bin. The value of the deviation corresponding to 68% of 
the area, about the nominal fit prediction, under the distribution is taken as the asymmetric 
statistical uncertainty on the fit. In the signal regions, the statistical uncertainty on the fit 
rises from 5% at the lower edge of the Tfx range to 17% at the limit of no data for Ajet > 3, 
and from 37% to 67% for Ajet > 8. 

A systematic uncertainty is assigned due to the range chosen for the control regions. 
To estimate the uncertainty on the nominal choice, the data are fit in all eight possible 77x 
control regions by increasing and decreasing the fit range by 0.1 TeV, and shifting it by 
±0.1 TeV. Each fit is used to predict the number of events in each 77x bin. The control 
region predicting the largest number of events and the control region predicting the smallest 
number of events for each i7x bin provide an estimate of the asymmetric uncertainty due 
to the choice of control region. In the signal regions, the systematic uncertainty due to the 
choice of control region rises from 2% at the lower edge of the 77x range to 8% at the limit 
of no data for Ajet > 3, and from 28% to 53% for Ajet > 8. 

To estimate the uncertainty in the analysis due to the choice of fit function, alternative 
fit functions are considered. Alternative fit functions are chosen such that in multi-jet 
simulated events (PYTHIA 8 and Herwig++) they provide a good fit in the control region 
and the extrapolation of the fit provides a good description of the simulated events in the 
corresponding signal region. In addition, the function should also provide a good description 
of the data in the control region. Only functional forms that fulfil these criteria for at least 
one inclusive jet multiplicity region are considered, and these are: 



(7.1) 
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(7.2) 

Po{l - xf^xP^^'^^ , 
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The systematic uncertainty due to the choice of fit function is assigned as the envelope of 
all alternative fit functions around the nominal function when fit to data in the control 
region. Figure 3 shows the alternative fit functions for different inclusive jet multiplicities. 
In the signal regions, the systematic uncertainty due to the choice of fit function rises from 
10% at the lower edge of the range to 46% at the limit of no data for > 3, and 
from 6% to 10% for Njet > 8. 

The jet energy scale (JES) uncertainties and jet energy resolution (JER) uncertainty 
are used in the MC validation methods and applied to the MC signal samples for the model- 
dependent limit calculations only. In the signal regions, the JES uncertainty rises from 1% 
at the lower edge of the range to 22% at the limit of no data, and the JER uncertainty 
rises from 0.6% to 2% for iVjet > 3. For A(jet > 8, the JES uncertainty rises from 21% at 
the lower edge of the range to 34% at the limit of no data, and the JER uncertainty 
rises from 7% to 11%. 


8 Results 

Figure 4 shows the extrapolation of the fits in the control region to the signal region with 
all uncertainties included. No data events are observed above Rt = 4.3 TeV, in agreement 
with the background estimate. 

To test the consistency of the data with the null-hypothesis (background-only hypo¬ 
thesis) a hyper-test statistic t = — ln[p-value™™] is defined where p-value™™ is the minimum 
local p-value in any inclusive A(jet and Rt region [42]. The Rt regions can range from a 
single 0.1 TeV bin to the entire range containing data. The hyper-test statistic takes ac¬ 
count of the trials factor (look-elsewhere effect). The most significant discrepancy in the 
observed signal region distributions is an excess in the interval 3.2 TeV to 3.9 TeV, for 
.Vjet > 4. This enhancement corresponds to a local p-value of 0.0043 which corresponds 
to a significance of 2.6 standard deviations [43] compared to the most probable value for 
the null-hypothesis of 2.3a. The corresponding f-value for the data is 5.4, and for the null- 
hypothesis the probability to find a value equal or greater than 5.4 is 0.4. This test shows 
that no significant excess is observed beyond the SM expectations for all choices of Rt and 
inclusive Vjet signal regions. 

Since black holes or string balls are likely to appear as an enhancement in the tail of Rt 
distributions, rather than as resonances, production upper limits are set in bins of inclusive 
H'j' (R™™), rather than Hj'. The predicted number of events in each 77^™ bin in the signal 
region is obtained by integrating the fit function from the Rt din of interest up to the 
kinematic limit of 8 TeV. The same integral is performed for the maximum and minimum 
number of predicted events obtained from the statistical uncertainty and each systematic 
uncertainty. The differences between the maximum (or minimum) number of predicted 


10 




1.5 2.0 2.5 3.0 3.5 4.0 4.5 5.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5 5.0 

H.^ [TeV] H.^ [TeV] 



1.5 2.0 2.5 3.0 3.5 4.0 4.5 5.0 


Ht [TeV] 



1.5 2.0 2.5 3.0 3.5 4.0 4.5 5.0 


Ht [TeV] 



Figure 3: The distributions, showing alternative ht functions for different inclusive 
jet multiplicities A^jet- The systematic uncertainty due to the choice of ht function versus 
exclusive H'y is given by the envelope of all the functions. 
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Figure 4: The H^' distributions, showing the data and extrapolated fits from the control 
region 1.5 < Ht' <2.9 TeV into the signal region H'j' >3.0 for each inclusive jet multiplicity 
Vjet- The uncertainty on the data points in both the distribution and ratio are due to the 
statistical uncertainty on the data only. The uncertainty band includes all uncertainties 
on the background prediction. Also shown are the expected black hole signals for three 
parameter sets of the CHARYBDIS2 non-rotating black hole model. 
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events and the nominal number of predicted events for each uncertainty are added linearly 
to obtain the total uncertainty on the number of predicted events in each bin. The 

uncertainties are added linearly since each is obtained from the same data in the control 
region. 

Using the observed number of events in each signal region compared to the predicted 
number of SM background events, model-independent limits on the observation of new 
phenomena at high and inclusive A^jet are set. In addition, model-dependent limits on 
several black hole and string ball signal models are set using the previous information, and 
estimates of the acceptance and efficiency for each model. The statistical procedure em¬ 
ployed uses the one-sided profile-likelihood test statistic [44]. The uncertainties are modelled 
with a convolution of Gaussian probability density functions describing the uncertainties 
on the signal or on the background. The upper limits are derived from pseudo-experiments. 

Counting experiments with as a function of are performed for each 

inclusive jet multiplicity Ajet > 3 to A^jet > 8. Using the number of data events, estimated 
background, estimated uncertainty on the background, luminosity, and uncertainty on the 
luminosity, upper limits are calculated on the number of events divided by the integrated 
luminosity. Upper limits on the visible cross section (number of events divided by the 
luminosity or cross section times acceptance times reconstruction efficiency) at the 95% CL 
are shown in figure 5. 

The expected uncertainty bands narrow at high Lfx as the test statistic becomes 
discretely distributed due to the extremely small background prediction. These model- 
independent limits on the cross section times acceptance times efficiency are as low as 
0.14fb at the 95% CL for minimum Lfx values above 4.3 TeV where no data events are 
observed. 

To set limits on various models, the visible cross sections are divided by the recon¬ 
struction efficiencies to obtain limits on fiducial cross sections defined at the particle level. 
The reconstruction efficiencies are calculated by taking the total signal efficiency times ac¬ 
ceptance and dividing by the fiducial acceptance. The fiducial acceptance at the generator 
particle leveP is defined from the simulated signal events with final states that pass the 
jet pt > 50 GeV and \ri\ < 2.8 requirements. The events are then counted in the different 
inclusive A(jet and Lfx regions. Detector resolution can cause migration of signal events to 
different regions of A(jet or iLx- This can cause the reconstruction efficiency to be larger 
than unity in some (A(jet, L^t) regions. 

The reconstruction efficiency depends on the particular black hole signal production 
model and kinematic region (A^jetj Lf™™) of interest. The efficiencies are determined over the 
range of n, Md, and Mth shown in figures 6-8. The efficiencies are practically independent 
of the number of extra dimensions in the model. The mean reconstruction efficiencies are 
about 88% with a variation between models of about 1%. The RMS spread in efficiencies 
for a particular model is about 4% with a variation between models of about 1%. 

The limit on the cross section times acceptance is obtained from the visible cross section 


^This includes parton showering and jet clustering, using the anti-fct algorithm with R = 0.4 on stable 
particles (cr > 10 mm), including electrons and photons, but excluding muons and neutrinos. 
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Figure 5: Upper limits on the visible cross section (cross section (rr) times acceptance (A) 
times efficiency (e)) at the 95% CL versus inclusive for different inclusive jet multipli¬ 
cities. The solid (dashed) lines correspond to the observed (expected) upper limits. The 
green (dark) and yellow (light) bands represent one and two standard deviations from the 
expected limits. 
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and the average reconstruction efficiency from many models, parameters, and kinematic 
properties. The limit on the cross section times acceptance is as low as 0.16 fb at the 
95% CL for minimum values above about 4.3 TeV. This is comparable to the results 
in ref. [17] in which missing transverse momentum is included in the definition of if it 
is greater than 50 GeV. The results presented here are also compared to another ATLAS 
analysis [14], in which a high-px lepton was required. The intersection of these limits with 
theoretical predictions for the cross section within the fiducial selections used in this analysis 
could be used to constrain other models of new physics resulting in energetic multi-jet final 
states. 

Exclusion contours are obtained in the plane of M/j and Mth for several models. For this 
purpose, counting experiments are performed to set a 95% CL cross-section upper limit for 
each signal model and each (Mth, M^i) value of that signal model in this analysis. In setting 
these limits, both the estimated background uncertainty and signal uncertainties are taken 
into account. The background and its uncertainty are the same as described previously for 
the model-independent limits. The uncertainties on the signal include the uncertainty due 
to the jet energy scale and jet energy resolution, the statistical uncertainty on the signal MC 
samples, and the uncertainty on the integrated luminosity. As the cross section for black 
hole production is known only approximately and is highly model dependent, no theoretical 
uncertainty on the signal cross section is applied. 

For each grid point in the Mth-M^i plane, the signal region which gives the lowest 
expected p-value is used. The most sensitive signal region for a particular signal model 
follows the kinematics of the signal model. For example, high-multiplicity signal regions 
are best for high-multiplicity signal samples, and high-77“™ regions are best for high-Mth 
signal samples. Observed and expected exclusion contours for different CHARYBDIS2 black 
hole models are shown in figure 6. The observed and expected exclusions for n = 2, n = 4, 
and n = 6 are shown. In each exclusion figure, lines of fixed Mth/Mj) (defined as k) are 
shown. The assumptions of the models are not valid for /c = 1, but are valid for /c ;§> 1. 
These lines therefore form useful guidelines as to the validity of the models across the plane. 
The results for non-rotating and rotating string balls are shown in figure 7. 

The exclusions tend to be stronger for higher n, due to the larger signal cross sections. 
For low values of Mth/M/r, where there are the fewest Hawking emissions (and where the 
semi-classical production assumptions are least valid), the limits worsen. The exclusions 
for non-rotating and rotating black hole models, with all other parameters identical, appear 
similar. Including initial-state radiation reduces the cross section and hence the (Mth, M/r) 
exclusion reach. A low-multiplicity remnant state weakens the exclusion reach for n = 2 at 
low values of Mth/M^), due to the reduced number of jets. For string balls, the exclusion 
for the rotating case is similar to the non-rotating case, in contrast to the result in ref. [14]. 

The exclusion contours for BlackMax models are shown in figure 8. They show the same 
general features as the ones obtained with samples generated by CHARYBDIS2. BlackMax 
uses a final-burst remnant model, which gives high-multiplicity remnant states [33]. Com¬ 
paring non-rotating and rotating CHARYBDIS2 and BlackMax results, shows this analysis is 
insensitive to the different remnant models. The results for the BlackMax model of produc¬ 
tion losses to photons is comparable to the results for the CHARYBDIS2 model of production 
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Figure 6: Exclusion contours in the plane for different black hole models in two, 

four, and six extra dimensions simulated with CHARYBDIS2. The solid (dashed) lines show 
the observed (expected) 95% CL limits. Masses below the corresponding lines are excluded. 
Lines of fixed (defined as k) are shown. The assumptions of the models are valid 

for k 1. 


losses to gravitons. Graviton emission in non-rotating black hole models weakens the ex¬ 
clusion slightly, as a greater number of decay products carry missing energy and do not 
contribute to the number of jets or iVx- 

Contour limits of Mth versus Mjy are presented for a variety of models. These limits 
can be interpreted in terms of lower-mass limits on black hole and string ball masses that 
range from 4.6 to 6.2 TeV. This is again comparable to the results in ref. [17] with the limits 
here being about 0.1 TeV higher in mass. The results presented here are also compared with 
those of ref. [14]. In the low-Mu region the results are comparable, while in the high-M^ 
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Figure 7 : Exclusion contours in the Mth-Mg plane for non-rotating and rotating string ball 
models simulated with CHARYBDIS2. The solid (dashed) lines show the observed (expected) 
95% CL limits. Masses below the corresponding lines are excluded. Lines of fixed Mth/Mg 
(defined as k) are shown. The assumptions of the models are valid for /c ^ 1. 


region the results presented here are a significant improvement over those in ref. [14]. The 
latter analysis is affected by a significant loss in sensitivity for the cases of rotating black 
holes and string balls, while the results presented here, and those in ref. [17], are rather 
independent of rotation. 
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Figure 8: Exclusion contours in the plane for different black hole models in two, 

four, and six extra dimensions simulated with BlackMax. The solid (dashed) lines show the 
observed (expected) 95% CL limits. Masses below the corresponding lines are excluded. 
Lines of fixed Mti^/Mo (defined as k) are shown. The assumptions of the models are valid 
for k 1. 

9 Conclusion 

The production of events with multiple high-transverse-momentum jets is measured using 
20.3 fb“^ of proton-proton collision data recorded at ^/s = 8 TeV with the ATLAS detector 
at the LHC. No significant excess beyond SM expectations is observed, and upper limits on 
the visible cross sections for non-SM production of these final states are set. Using models 
for black hole and string ball production and decay, exclusion contours are determined as 
a function of mass threshold and the fundamental Planck scale. 
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The limit on the cross section times acceptance can be obtained from the visible cross 
section and the average reconstruction efficiency taken over many models, parameters, and 
kinematics. The limit on the cross section times acceptance is as low as 0.16 fb at the 
95% CL for minimum values above about 4.3 TeV. 

Contour limits of Mth versus Mjy are presented for a variety of models. These limits 
can be interpreted in terms of lower-mass limits on black hole and string ball masses that 
range from 4.6 to 6.2 TeV. 
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